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Abstract. Modified Newtonian Dynamics is an empirical modification to 
Poisson's equation which has had success in accounting for the 'gravitational field' 
4> in a variety of astrophysical systems. The field <E> may be interpreted in terms of 
the weak field limit of a variety of spacetime geometries. Here we consider three 
of these geometries in a more comprehensive manner and look at the effect on 
timclike and null geodesies. In particular we consider the Aquadratic Lagrangian 
(AQUAL) theory, Tensor- Vector-Scalar (TeVeS) theory and Generalized Einstein- 
Mther (GEA) theory. 

We uncover a number of novel features, some of which are specific to the theory 
considered while others are generic. In the case of AQUAL and TeVeS theories, 
the spacetime exhibits an excess (AQUAL) or deficit (TeVeS) solid angle akin to 
the case of a Barriola-Vilenkin global monopole. In the case of GEA, a disformal 
symmetry of the action emerges in the limit of V$ — ► 0. Finally, in all theories 
studied, massive particles can never reach spatial infinity while photons can do so 
only after experiencing infinite redshift. 



1. Introduction 

Einstein's theory of General Relativity is a supremely successful theory on scales of the 
solar system and below. Although specific predictions about such diverse phenomena 
as the gravitational redshift of light, energy loss from binary pulsars, the rate of 
precession of the perihelia of closed orbits, and light deflection by the sun are not 
unique to General Relativity, it must be regarded as highly significant that General 
Relativity is consistent with each of these tests and more. 

The success of General Relativity on larger scales still is perhaps less clear. On 
the largest observable scales, observations appear to have converged upon a simple 
model: the so-called Concordance Model. The addition of a cosmological constant 
and an approximately pressureless fluid (termed Cold Dark Matter) to a system of 
General Relativity coupled to the fields of the standard model of particle physics 
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(and in conjunction with appropriate initial conditions) appears sufficient to give an 
excellent account of observables such as the large scale structure in the universe and 
the cosmic microwave background anisotropies. 

Despite the success of the Concordance Model, the requirement of the existence of 
an extra, stable matter field as yet unobserved in any small scale experiment and the 
accompanying requirement of a (by all conventional standards) very small but non- 
zero cosmological constant has led some to question whether there is something lacking 
about the model. A differing model of the universe at the level of linear cosmological 
perturbations would presumably amount to a different account at the level of nonlinear 
perturbations: i.e. on small scales such as galaxies or the solar system. Crucially 
though, galactic scales do indeed exhibit mass discrepancies: stars move as if pulled 
on than a mass greater than the observed luminous mass, gravitational lensing by 
many galaxies is greater than one would expect assuming General Relativity and 
the observed luminous mass. Currently it is not entirely clear that they are of the 
precise nature predicted by the Concordance Model. In light of the uncertainty about 
what constitutes the 'dark sector' in the universe, and the possible limitations of the 
Concordance Model, one may instead attempt a phenomenological description of the 
mass discrepancies on small scales and attempt to work from there to an appropriate 
cosmological model. 

Indeed, in 1983 it was realized by Milgrom [I] that the emerging evidence for the 
presence of dark matter in galaxies could rather follow from a modification either to 
how 'baryonic' matter Q responded to the Newtonian gravitational field it created or 
to how the gravitational field was related to the baryonic matter density. Collectively 
these ideas are referred to as MOdified Newtonian Dynamics (henceforth MOND). By 
way of illustration, MOND may be considered as a modification to the nonrelativistic 
Poisson equation [2] : 



where $ is the gravitational potential, V is the grad operator, ao is a number with 
dimensions Length -1 , Gjy is Newton's constant and p is the baryonic matter density. 
The number a is determined by looking at the dynamics of visible matter in galaxies 
and a value of the order of the inverse Hubble length seems to be preferred [3] . The 
function fJ,(x) would simply be equal to unity in Newtonian gravity. In MOND, the 
functional form is only fixed at its limits: /i — > 1 as x — > oo and /i — > x as x — ¥ 0. We 
will call the limit /i — > x the 'modified gravity' (MG) regime. 

The success of ([1]) may simply reflect the behaviour of dark matter on these scales. 
Alternatively, the complete absence of anything other than a gravitational potential 
and baryonic matter from this equation may appear to indicate that the origin of (JXJ) 
lies along the lines of a modification to how the baryonic matter affects spacetime i.e. 
a modification to gravity. 

Considerable work has gone into exploring theoretical implications ( [H [7J [5J 
M HDD and observational implications ([TTl[Illl[Il[Il[I5J[Il[n])of equation 
If MOND is taken to be a modification to gravity, it would seem necessary to find a 
theory where ([IJ emerges from a relativistic framework, or a wider framework wherein 
at least the successful predictions of General Relativity could be recovered in the 




(1) 



X Baryonic matter is denned to be matter which consists of ingredients from the standard model of 
particle physics 
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appropriate limit. In this absence of this, it is unclear precisely when one is able to use 
MOND and when it is appropriate to compare the predictions of MOND with those 
of the dark matter paradigm. This was quickly recognized and in 1984 Bekenstein 
and Milgrom introduced a generally covariant theory which could reproduce MOND 
phenomenology in the absence of dark matter [2] . We refer to this theory as AQUAL 
(not to be confused with a nonrelativistic Lagrangian theory sometimes given the same 
name), named after the AQuadratic Langragian of the theory. 

Since then, a variety of other relativistic models have been proposed which seek to 
reproduce aspects of equation ^ in the quasistatic, weak field limit, (see for instance 
(TBI HSl 1201 EH H3 HSl H31 HSl HSl I2Zl Of these, we single out two theories 

which we will further use in this work, namely Bekenstein's Tensor- Vector-Scalar 
(TeVeS) theory [2U] and the generalized Einstein-iEther(GEA) theory [35]. A feature 
common to all these models is that the transition from Newtonian gravity to modified 
gravity is accomplished via the introduction of a rather unwieldy function into the 
model's action which is perhaps an indication that a different approach is called for. 
Furthermore, none of the models can claim to be as successful as the Concordance 
Model on cosmological scales [29l [3Q1 [31] . 

In this paper we adopt a more conservative viewpoint: It may be that if there is 
underlying 'new physics' to the MOND paradigm, it may not be convenient to express 
its equations in the form (p}. We will simply assume that the new physics allows for 
a modified gravity regime which for appropriate sources is described at the largest 
distances from the source by a static, spherically symmetric spacetime metric. We are 
interested in the properties of this asymptotic spacetime. For instance, this spacetime 
may be expected to approximate that due to a black hole or a massive bound object at 
sufficiently large distances. This approach will enable us to deduce general properties 
of these spacetimes and need not rely on the precise details of how the transition from 
General Relativity to modified gravity is accomplished or the details of the underlying 
theory. 

2. The three relativistic MOND theories 

We consider three distinct relativistic MOND theories. For each theory, we will provide 
the action and a brief description regarding the form of the free function involved, 
followed by the relation of its parameters to the non-relativistic limit. We will then 
provide the spherically symmetric, vacuum field equations and proceed to solve them 
in the limit of large distance from the source. This is the limit where the 'modified 
gravity' regime is appropriate. 

In all cases we adopt a static and spherically symmetric metric. In most of the 
cases we employ isotropic coordinates {t, r, 9, (p} for which the metric takes the form 

ds 2 = _ e 2* dt 2 + e 2* dL 2 (2) 

where dl? = ^ijdx % dx^ — dr 2 + r 2 {d9 2 + sin 2 9d(p 2 ^ is a flat metric in spherical 
coordinates and $(r) and ^(r) are two functions of r which are to be determined by 
the field equations for each theory. We note that both $ and \P are determined up-to 
a constant, i.e. the transformation $ — > $ + $o and — > + for constants $o and 
"Jo can be undone by a further coordinate rescaling of t and r. As we show further 
below, the metric for the three theories in the modified gravity limit, takes the form 

/ \ 2c / \ 2c 

AQUAL: ds\ = - ( — J dt 2 + dl? (3) 
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TeVeS : ds 2 T 
GEA : ds% 

for a dimensionless number c and a scale ro emerging as an integration constant which 
can be chosen to roughly mark the onset of the modified gravity regime. 

The use of isotropic coordinates makes it straightforward to determine the non- 
relativistic limit by simply expanding e 2 * — > 1 + 2$ + . . . and e 2 * —¥ 1 + 2^ + . . .. 
For all the three cases above, we find that $ = c In j- which is precisely the MOND 
potential. This enables us to identify the dimensionless number c with ^/Uj^Moq (M 
being the mass of the source). To get a feeling for the order of magnitude of this 

number we can write it as c ~ 1.4 x 10~ 12 y^^, where Mq is the solar mass [3]. This 

is a tiny number by all standards, except when M starts to approach the total mass 
within our Hubble volume. The non-relativistic limit is thus expected to hold for most 
cases of interest. It is still desirable however to look at the full relativistic metric. For 
instance, while 4> is the same for all three theories (as required for reproducing the 
MOND limit), the potential * can be different. In particular AQUAL gives * = $ 
while TeVeS and GEA give *S? = — This is a result that could not have arisen 
from ([T]). As we shall see further below, novel features can arise by the use of the full 
relativistic solution. 



ro 



- dt 2 + - 



In 



To 

i 6 



dL 2 



dt 2 



dL 2 



■In 



(4) 
(5) 



2.1. Aquadratic Lagrangian (AQUAL) gravity 



The simplest of the theories we consider is the AQUAL theory of Bekenstein and 
Milgrom [2]. The theory depends on an Einstein metric g a b and a scalar field <f>. The 
physical metric g a t, is given by a conformal transformation g a t, — e 2 ^g a b- The action 
is 



S A [g a ,<t>] 



l 



16ttG 



R-TsfW + ^matter [<?] (6) 



where £q is a scale and where 

0% ~abx 



X = e 2 g ab V a cbV b 



(7) 

The theory depends on a free function f(X) which must have appropriate limits to 
recover Newtonian or MOND gravity in the non-relativistic case. In particular, the 
non-relativistic limit of the theory gives back the MOND equation (TTJ) with /x = 2 +f x 



where f x = Jf. Defining limx-s-oo f(X) = f X and lim x ^ /(A) 
find that the Newtonian limit fixes the Newton's constant Gn 
4> c is the background (cosmological) value of and gives the Parameterized Post- 
Newtonian(PPN) parameter 7 = In the MOND limit, we find that Milgrom's 



/1A 3 / 2 

= 2±£o e 20 cG where 
Jo 



constant ag is given by ap — 



4/0 



_ 3/i(2+/ )^ to 4 

We find the field equations assuming isotropic coordinates. The Einstein metric 



is given by 500 = — e 2 * with $ = <f> + cj) and (jij 
scalar field equation for our metric ansatz is 

d 



2* 



dr 



$+3' 2 

3 V 



]x Tr 



= 



7ij with ^ = 4 r + 0. The vacuum 



(8) 
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where X = t\e 2 * ' 4>' 2 , while the Einstein equations are 

2*" + *'(*' + J)=-^/e 2 * (9) 
H ($' + + (2$' + §') = ie 2 * (x/x - -/) (10) 

$" + §" + I + + ($')2 = -^l/e 2 * (11) 

In the MG limit / = /iA 3//2 so the scalar field equation integrates once to 



2C e 2 ce 



(12) 



where Co is an integration constant, related to the mass of the source. As r — > 00 the 
RHS of the Einstein equations goes as X 3 / 2 e 2 * — > e - 3 */ 2 -* which will decay faster 
than any of the terms in the LHS. This means that in the MG limit as r —> 00, 
the Einstein-frame metric must be a solution to the vacuum Einstein equations 
and for a static spherically symmetric spacetime that solution is unique; it is the 
Schwarzschild solution. We may thus write $ = ^ + ... and = — + 

while <fi' = ~ (l — r£ + . . .) where c = 3 2 ^ p . The last relation integrates to 

4> = cln ~ + + . . ., where the integration constant is ignored as it is irrelevant. 
This gives the physical metric as 

2 c 



r 



ds z = I — 

' - (2 7 )$1 +-)^ a ] (13) 

which is the metric <J^5j) in the limit r — > 00. 

A comment is in order. In constructing (1121) we have retained the positive sign 
when taking the square root. The reason is that c < is unphysical. To see this, 
we take the non-relativistic limit of ([T3|) which gives $ = —\c\ lrx(r/ro). This implies 
that the acceleration a = — V$ = |c|- points away from the centre which would 
render gravity being repulsive (this is the analogue of the negative-mass Schwarzschild 
solution). 

It is instructive to find the physical metric in Schwarzschild coordinates f . The co- 
ordinate transformation in the r — > 00 limit implies that = ( ^- J (l 

which can be inverted to give ±- = ( £-) 1+ " \l + „, 1 c f\ (Q)^ + .. 

b r V'' / L 2(l+c)r \ r/ 

f — > 00 as r — > 00. Thus we find the metric in Schwarzschild coordinates as 



1r ~ •••) 
therefore 



ds z = - [ — ) ^ dt 2 + -^T7 + ^ 2 ( 14 ) 
\r ) (1 + c) 2 

The factor appearing in the df part of the metric is quite important. It implies 
an excess solid angle given by 1 + c. In fact, the spatial metric is of exactly the same 
form as the Barriola-Vilenkin global monopole solution [32] (the original solution has 
a deficit angle) so this MONDian spacetime can be thought as a warped Barriola- 
Vilenkin global monopole. 
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The divergence of the conformal factor appearing in the metric may be a source 
of concern. However, no curvature singularity arises. The curvature tensors go as 

R ab cd = pa& crfr -2(l+ C ) ^ _ pa r 2(l+c) R = _g c ^ + ^ r _ 2 (l+c) ) where the p_ 

tensors are constant. So forming invariants by contracting various types of curvature 
tensors together, or their derivatives, we would get that a typical invariant would be 
I n _d = 2Ti(i+ C )+d ■> where J is a real number, n is the number of curvature tensors 
involved and d is the number of derivatives acting on the curvature tensors. Thus all 
curvature invariants vanish as r — > oo. 



2.2. Tensor-Vector-Scalar (TeVeS) gravity 

Bekenstein's TeVeS theory combines the AQUAL theory with a unit-timelike vector 
field A a considered by Sanders 33J in the MOND framework on phenomenological 
grounds. The relation between the Einstein and physical metric is changed to 
gab = e~ 2 ^g a b — sinh(20)^4 a A(,, the so-called disformal transformation. We shall adopt 
the diagonal- frame formulation of the theory (see appendix of [27]) where the scalar 
field is decoupled from the vector field. The action takes the form 

K abcd V a A b V c A d 

(15) 



S T [g, <p, A] = S A [g, <j>\ - / <^xy/^§ 



X(A a A a + 1) 



where A a = g ab Ab and the tensor K abcd is given by 

K abcd = Cig ac~bd + C2 gab~cd + ^ad-bc + ^a^-bd^ ^g) 

for a set of constants {C1...C4}. This formulation of the theory is more general 
(corresponding to the proposal in [34]) than the one given by Bekenstein which 
corresponds to c\ = 2K — |, c 2 = —5, c 3 = — 2K + | and C4 = K — j, where K is 
Bekenstein's constant in [20] . As it happens, spherically symmetric solutions depend 
only on the combination c\ — C4 and not on c 2 or C3. Without loss of generality we 
shall therefore set ci — C4 = K from now on. 

Once again, the theory depends on a free function f(X) which must have 
appropriate limits to recover Newtonian or MOND gravity in the non-relativistic 
case. In particular, the non-relativistic limit of the theory gives back the MOND 
equation |T]) with /i = ( 1 + ^j^-J 2 -k+f x ' Defining Yimx^oo f(X) = foX and 
limx-i-o /(A) = /1A 3 / 2 we find that the Newtonian limit fixes the Newton's constant 
G N = 2=k) G and § ives the PPN P arame ter 7=1. In the MOND limit, we 

find that Milgrom's constant ao is given by ao — S~ 1 r t— ^ c t-. 

We find the field equations assuming isotropic coordinates. We assume that the 
vector field has no spatial component. The Einstein metric is given by goo = — e 2 * 
and (jij — e jij. The unit-timelike constraint A a A a = —1 then fully determines the 
vector field as A a = (e*, 0) and A a = (-e~*, 0). We further find that $ = $ + and 
* = <l - 4>. 

With this in hand, the vector field equations arc identically satisfied. The vacuum 
scalar field equation for our metric ansatz is unchanged from the AQUAL case, while 
the Einstein equations become 



2*" + #' ( f + = -\fe™ - K 



r) 111* 



r , 



(17) 
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2 / , 



+§' (2l>' + \E' / 



' 



$" _|_ i|r" . 



2* 
1 

'2? 



Xf 



x 



1 -/ 
2 



1 



12 



(18) 
(19) 



In the MG limit we have that / = f\X 3 / 2 , hence the scalar field equation can be 
integrated as in the AQUAL case. Once again, in the r — > oo limit, the scalar- 
dependent terms in the Einstein equations vanish more rapidly than the metric terms 
and the corresponding field equations become those of Einstein- ./Ether theory. The 
solutions in that case have been found by Eling and Jacobson [35]. We do not need 
the exact solution here but only the fact that as Eling and Jacobson show, the solution 



is asymptotically flat. We may then expand <f> 



cln ■ 



2r 



as 



. where c > is given as before. 
(2 



- . . . and * = + while 

r 1 

This gives the physical metric 



ds 2 = - — 







r 


r 


ro 






ro) 







1 



1 



C )$! 



(2 + c)$i 



dt 2 



dh l 



(20) 



which is the metric (j4]) in the limit r — > oo. 

It is instructive once again to find the physical metric in Schwarzschild 
coordinates f. Assuming for the moment that c < 1 (we shall motivate 
this further below), the coordinate transformation in the r — > oo limit implies 

(2+c)3>i , 



that -f = 

^0 



(ft) 



1 - 



which can be inverted to g: 



lve — 

ro 



\r 



I , (2+c)$ 



2(l-c)r V f 



therefore 



the metric in Schwarzschild coordinates as 



ds 2 



ro 



dt 2 



df 2 



r — > oo as r — > oo. 



- 2 dn 2 



Thus we find 



(21) 



(1-c) 2 

As in the case of AQUAL, the spatial metric is of the same form as the Barriola- 
Vilcnkin global monopole solution [32] , the factor appearing in the df part of the 
metric implying now a deficit solid angle given by 1 — c. 

It is interesting that as c — > 1 the transformation becomes singular. Does it 
even make sense if c > 1? We can answer this by first looking at the curvature 



The curvature tensors go as R 



nb 



cd 



P 



ab „2(c-l) 
cd 1 ' 



invariants. 

R = 2(1 — c)r 2 ^ c-1 \ where the P-tensors are constant. So forming invariants by 
contracting various types of curvature tensors together, or their derivatives, we would 
get that a typical invariant would be = r2 „ (1 '{ c)+d , where J is a real number, n is 
the number of curvature tensors involved and d is the number of derivatives acting on 
the curvature tensors. Thus all curvature invariants vanish as r — > oo provided c < 1 
and if c > 1 they diverge. This indicates that c > 1 is not physical. Another way to 

see this is as follows. The proper distance is D — (ro) ' nence ' for c > 1 we 
find that D — > as r — > oo. Therefore, if the solution was valid for c > 1 there would 
be a naked singularity at zero proper distance, surrounded by MONDian spacetime 
up-to some finite distance ro . We conclude that the case c > 1 is not physical in TeVeS 
theory. 



The Geometry Of Modified Newtonian Dynamics 



8 



2.3. Generalized Einstein Mther gravity 



The Generalized Einstein Aether theory involves a unit-timelike vector field much like 
TeVeS, however, it has no scalar field. The theory is formulated with a single metric, 
and the MOND behaviour is recovered by using a noncanonical kinetic term [22] for 
the vector field. Its action is as follows: 



S[g a \A a ,X] 



1 



16ttG 



d 4 x. 



-fJ 



R 



+ A (A a A a + 1) 



^matter id] 



where 



X = £' z K abcd V a A b V c A d 



(22) 



(23) 



and K is given by (1161) with g —> g. 

The theory depends on a free function f{X) which must have appropriate limits to 
recover Newtonian or MOND gravity in the non-relativistic case. The non-relativistic 
limit of the theory gives back the MOND equation (JTJ) with /i = 2 ^_^ x , where 
K = C\ — C4. The Newtonian limit is recovered when limx-s-oo f(X) = X in which 
case the observed Newton's constant is fixed as Gn — t^k^ ano - * ne PPN parameter 
7=1. To get MOND behaviour we define a new function h(X) = f(X) — so that 
fi = MOND is recovered if lim.x->oh(X) = fi(-X) 3 / 2 in which case we get 

that Milgrom's constant ao is given by oq — ^ . 

We now proceed to find the relativistic MOND solution. We assume that there 
is no spatial tilt to the vector field A a . Consequently, the vector field is entirely fixed 
in terms of the metric by the fixed-norm constraint and the vector field equations 
are identically satisfied. After calculation it may be shown that the vacuum Einstein 
equations reduce to the following: 



(hx&)' + h x & 



r 



(Xh x -h)=0 



and 



2 ($" + *") + -($' + + + * 



'\2 



.,2* 



Xh 



x 



(24) 



(25) 



Substituting h = fi(-X) 3 / 2 for pure MOND, where X = 
become: 



1 







and 



where Y = $ ■ 



-Y' + -Y' 2 = 
r 2 



Y" 

The solution is 
$ =$o + 3 In (l 

Y = Y n + 2 In I 1 



c r 
-In — 

3 r 



Iro 
2 r 



ATe the equations 
(26) 

(27) 

(28) 
(29) 
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where we have the constants of integration c,$o? r 0i an d Y . As can be seen in section 
(|3.2p , the number c has a clear physical interpretation. The numbers $ an d Y m ay 
without loss of generality be set to zero by rescale the temporal and radial coordinates. 
We this arrive at the metric 



ds 2 £ = 



1 + - In — 



3 \r 



6 h — lo.\ 2 

L 2r ' 



dt z + i *U. dL 2 (30) 



1 + # In X. 

3 V ro 



For large distances r > ^ we recover the metric ([S]). 

The curvature tensors go as R ab cd = ^P ab cd {\n^) R\ = ^P\^^) R = 
^■P(ln^), where the P-tensors are polynomial functions of m^-. So forming 
invariants by contracting various types of curvature tensors together, or their 
derivatives we would get that a typical invariant would be I n .d = r 2^+d ■/(hi where 
J is a polynomial function of In n is the number of curvature tensors involved and 
d is the number of derivatives acting on the curvature tensors. Thus all curvature 
invariants vanish as r — > oo. 

Finally we note some novel behaviour at the point X — > 0. Here the Einstein 
equation (I24|) is identically zero leaving only equation equation (|25[) which becomes 
again a differential equation for the variable Y = $ + ^. It may also be shown 
that the action (|22p in this limit is only a function of Y and is thus invariant under 
$ — > <f> + j3(r) and $ — >• $ — This disformal symmetry can be more explicitely 

written as g a b — > e _2 ^ffa6 — sinh(2 f3)A a Ab and A a — > e~^A a where /3(r) is some 
function of r. Interestingly the metric transformation is of the same form as the 
transformation between Einstein and physical metrics in the TeVeS model [20) . 

3. Timelike Geodesies 

3.1. General motion of test particles on planar orbits 

We now seek to look at effects which do not rely on the nonrelativistic conditions 
being satisfied. 

Consider a test particle's four velocity U a . In terms of the proper time r along 
the trajectory we have U a = dx a /dr = x a . Furthermore, without loss of generality 
we may consider trajectories lying in 6 = tt/2. Hence, the unit-timelike condition on 
the particle's four velocity may be expressed as follows: 

-l = -e 2 *i 2 + e 2 *(r 2 +rV) (31) 

Furthermore we have some constants of the motion. The first is the projection of the 
four velocity along the timelike Killing vector (dt) a : 

E = -g ab (d t ) a U b = e 2 *t (32) 

The second is the projection of the four velocity along the angular Killing vector {d v ) a : 

L = 9a b (d v ) a U b = e 2 *r 2 ^ (33) 

These expression may then be put in the unit-timelike condition U a U a g a b = — 1 
obeyed by the test particle, which after rearranging becomes 

= i e 2$+2* -2 + e — + i_ e -2*+2* (34) 

We now define the new variable £ = J e* + *<ir, in terms of which the energy 
equation becomes: 
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E 2 

£ = — = T+V (35) 



2 

T ee ^ 2 (36) 



e 2 * 



Hence, the problem of motion of test particles is analogous to a one dimensional 
classical mechanics problem of a particle with energy per unit mass E 2 /2, kinetic 
energy per unit mass T, and potential energy per unit mass V(£). 

3.2. Circular Orbits 

In keeping with our classical mechanics analogy, circular orbits will exist given 
appropriate initial data and the existence of a value of £ which satisfies dV/dl; = 0. 
Assuming the existence of such a point, we may relate the angular momentum to the 
geometry at that point, leading to the following relation: 

d± r 3 

L 2 = e 2* d£ (3g) 

dr i „ I d<S> x 



We can express this in terms of a more familiar quantity, the norm of the spatial part 
of the particle's four velocity, defined as v 2 = h a f ) U a U b where h a b is the spatial metric. 
It may be checked that L 2 = v 2 r 2 e 2 * and so 

d± r 

•• 2 - (39) 



dr , 



( d£ _ d±\ 
\d£ dt ) 



Clearly then if e 2 * and e 2 * arc monomial functions of r then v will asymptote to 
a constant. This a fully relativistic realization of the so-called Tully-Fisher relation 
that MOND reproduces. In AQUAL we have v 2 = c = ^GMa . In TeVeS we have 
v 2 = c/(l — 2c). As we shall see, this result is to be expected: for c > (1/2) there 
exists no minimum of the potential V(£) and hence there cannot exist any circular 
orbits. Finally, we see that although the GEA model reproduces the TeVeS result in 
the nonrelativistic limit, this eventually becomes a poorer approximation and v 2 is no 
longer approximately constant across increasing radii of orbits. This transition will 
happen when ln(£/r ) becomes of order unity, representing a significant deviation of 
$' and fy' from power law behaviour. By inspection this will happen at a radius we 
will call r a which is found to be: 

r a = r e- (40) 

This is a novel feature. The AQUAL and TeVeS metrics maintain 
characteristically 'MOND' behaviour to infinitely large distances from the source, 
as is evidenced by v 2 remaining constant as one varies orbital radius. This is not 
the case however in the GEA model, which approximates TeVeS only in the regime 
ro <C r <C r a . The extent of this regime depends on the value of c = \fGMa~Q and 
thus the mass of the sourcing object. By way of example, the value of c exterior to 
a galaxy such as the Milky Way is of order 10~ 6 , meaning that here r a is a vastly 
greater value than the radius at onset of the modified gravity regime. 
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3.3. The Precession of almost-circular Orbits 

We now seek to generalize to examine the properties of orbits that are not circular. 
Differentiating the energy equation we have that: 

« = -f <«» 

In order to deduce how the orbital radius varies with angle, we would next like 
to convert £ into d 2 £/dp 2 . From the conservation of angular momentum we have 

7 J — 2'J' 7 

-f- = — -f-, therefore, 

dr r 2 dip ' 1 

(42) 

r 2 d<p \ r 2 dp J 

For simplicity of calculation we further define a new variable u — — j (e~ 2 * / 'r 2 )d£ and 
in terms of u rather than £, the orbital equation becomes 

d 2 u 1 dV 

d^=-J 2 ^ (43) 
A strategy for showing whether orbits are closed is to look at almost circular 
orbits, i.e. trajectories with u — u p + uo where uq is the function evaluated at the 
radius of a circular orbit and u p is small in some suitable sense. In this situation we 
may look to Taylor expand the right hand side of equation (1431) in powers of u p . The 
assumed smallness of u p means that we can, to a good approximation, neglect terms 
in u 2 and higher orders. The result is an equation where the second derivatives of 
u p with respect to p are proportional to u p . The solutions to this equation thus be 
Up = u P Q cos(i7(p + b) where u p q and b are constants of integration. If fi is an integer 
then almost circular orbits are closed. 

3.3.1. General Relativity For comparative purposes, we first demonstrate this 
calculation for the Schwarzschild solution of General Relativity. In isotropic co- 
ordinates the metric takes the following approximate form: 

ds GR = - [1 - + 2 *w + ■ • •] dt 2 + [1 + 2$ N + . . .} dL 2 (44) 

where $jv = QM., One may proceed with the method described above to obtain 
a solution of the form u p — u p oCOs(flGRP + 6), where there is a one-to-one 
correspondence between the field u p (r) and r. The quantity Qgr is found to be 

equal to ^Jl — rj 1 — d,G r M , where r c is the radius of the circular orbit that is 

perturbed around. Therefore, u p has an angular period of approximately 2n(l + 3< ^ ) 
radians. Significantly, the deviation from closedness decreases as the corresponding 
circular orbital radius increases. 

3.3.2. AQUAL Recall that the asymptotic metric in AQUAL is given by iJHJ. 
Therefore, £ is not identical to r but is rather £ = J e 2 *dr. However, there is a 
compensating simplification in the relation of the variable u to the original coordinates: 

-2* f 2*-2* -j 



It is then simple to write V as a function of u: 
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The solution for u is found to be: 



-(^) (47) 



r \crg 

After some straightforward algebra one finds that: 

^ = -2(l + cK (48) 

We see then that the angular frequency of u p in AQUAL is given by Vt 2 A = 2(1 +c). 
Therefore, near circular orbits in AQUAL are not closed unless the number c is an 
integer. Note also that the angular period is significantly different from 2ir, even 
if c <C 1 and in the general case of non-integer c, the deviation from closedness is 
independent of orbital radius. 

3.3.3. TeVeS We recall that at far distances from the source, the metric in TeVeS 
is given by (j3|. As in General Relativity, £ may be simply identified with r. We first 
find the variable u. It is found to be defined via the equation: 

Meanwhile the potential V may be expressed in terms of u simply as: 

V= i[(l-2c)r u]^r + ^_ (1 _ 2c ) 2 U 2 (50) 

As before we find uq, the value that u takes at the circular orbit. This is found to be 
as follows: 



1 l_ / (l~2c) L 2 

(1 - 2c) V Q \ cT 2 " 



After a straightforward but lengthy calculation it may be shown that the perturbation 
Up obeys the following equation: 

d 2 u 

f = -2(1 - 2c)(l - c)u p (52) 



dip 2 

Therefore the frequency is given by Qj, = 2(1 — 2c) (1 — c), which will typically be 
markedly different from integer values. As in AQUAL, the angular frequency differs 
significantly from 2ir even if c <C 1. It is noteworthy that for the range 1/2 < c < 1, 
the perturbation u p will grow exponentially. Indeed, by inspection the potential U(£) 
can no longer have a minima in this case, and will inevitably push the test particle to 
smaller values of £. Recalling that c = \/GMao, we see that in TeVeS if there exists an 
object whose mass is greater than l/(4aoG) and has a region exterior to it in the deep 
MOND regime, then this region should be expected to be devoid of matter. Similar 
conclusions will hold in the GEA model as long as the condition r <C r a is satisfied. 

3.4- Radial Trajectories 

We now ask whether it is possible for a test particle to be projected to spatial infinity. 
The method to do this with the least energy would clearly be to do so with vanishing 
angular momentum L as this will only detract from how much of the energy £ can go 
into the radial motion of the particle. 
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3.4-1. General Relativity Recalling the metric (1441) we may immediately write down 
the zero angular momentum energy equation: 

1 •„ If 2GM \ , . 

Consider a test particle propelled radially outwards from position £o and with 
velocity (£)o- For the trajectory to reach £ = oo having slowed down to £ = 0, it must 
have energy £ = 1/2, and of course is the least amount of energy necessary to reach 
spatial infinity. Additionally this implies that we have the relation (£)o = \/2GM/£o 

3.4-2. AQUAL In the AQUAL case, although £ is no longer identical to r, it may 
readily be checked that £ = oo indeed corresponds to spatial infinity for all positive 
values of c. The zero angular momentum energy equation takes the form: 



2 2 \r 0/ 

We immediately see that there exists no finite value of £ for which the particle may 
begin at a finite value £o and reach £ = oo, or equivalently, there is no finite initial 
velocity with which the particle could be propelled to each spatial infinity. 

3.4.3. TeVeS Recall that £ is identically equal to r. The proper distance D between 
£ = £0 an d £ = 00 is then: 

Since c < 1 in TeVeS, these co-ordinates span an infinite proper distance. The zero 
angular momentum energy equation is then: 



Precisely as in AQUAL, no finite £ permits the test particle to reach spatial infinity. 

3-4-4- GEA As in the case of TeVeS, £ may be identified with r. The proper distance 
between D between £ = £0 and £ = 00 is of the form: 

D= f°° (i + l]n(i.)) (57) 



?0 



3 \r 



This integral diverges for all positive values of c and so isotropic coordinates are again 
appropriate. 



S =2^ +H 1+ 3 ln UJJ (58) 
Once again, no finite £ permits the test particle to reach spatial infinity. Therefore, 
for all of the spacetimes considered here, there can exist no unbound orbits. Put more 
informally, in AQUAL, TeVeS, and GEA the escape velocity for all massive particles 
is the speed of light. 
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4. Null Geodesies 

We may assume that the trajectories of test photons follow null geodesies i.e. their 
wave vector k a = dx a /d\ (where A is an affine parameter) satisfies k a k a = 0. As 
before, we may make use of constants of the motion E and L which, as the notation 
suggests, correspond to constants of the motion associated with the projection of 
k a along the time like and d v directed Killing vectors respectively. The null norm 
condition is then simply: 

E 2 = e 2*+2^2 + ^e 2 *- 2 * (59) 

Therefore again we may make the analogy with a one dimensional classical mechanics 
problem of a test particle with a trajectory £(A) = J e* + *dr. Its energy equation is 
then £ = T + V with 

F 2 

(60) 

T = \i 2 (61) 

V = - L2 e 2$-2tt ( q 2 \ 

Jf.,1. Gravitational Redshift 

We now consider the effect on a photon's frequency as it passes through the spacetime. 
We consider two stationary observers, one at coordinate r\ and another at r-i where 
r2 > T\. Their four velocity U a (r) may be related to the spacetime's timelike Killing 
vector £° = (1, 0, 0, 0) as follows: 

U a = , 1 C = -^=C (63) 
V-9ab£, a £, b V-900 

Consider a photon trajectory with wavevector k a emitted at r\ and absorbed at 
ri . At T\ , the frequency of the photon measured by the stationary observed is simply 
uji = —k a U a (ri). Similarly, the frequency of the photon measured later by the observer 
at r2 is given by ui2 = —k a U a (r2). We may express the ratio of these frequencies in 
terms of the wavevector an the Killing vector using equation (|63p . yielding: 



^2 _ y/-goo(ri) fc a £ a (r 2 ) ^ 
w i ^-500(^2) k a ia{ri) 

However, from the geodesic equation and Killing's equation, it follows that 
C k b(^ a k a ) = and so the ratio ^"(^j is simply equal to unity. Thus the ratio 
of frequencies depends only on the varying norm of the Killing vector. 

4.1.1. TeVeS and AQUAL As the time-time component of the metric is identical 
in TeVeS and AQUAL, they will describe identical gravitational redshift. The result 
may simply be read off to be: 



W2 \ / n 



(65) 
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4- 1.2. GEA Again we may read off the result given the metric: 

/ I -L _ h, I LL \ \ " 

LU 2 



^1/ GEA 




(66) 



Clearly in all three models, u^/^i — > as r 2 — > oo. Therefore, no photon can 
reach infinity with non-zero energy. This an equivalent result to the case of massive 
particles. 

4-2. The Deflection Of Light 

Although the frequency of a photon diminishes to zero as it approaches infinite distance 
from the source, null geodesies can be extended to these distances for any positive value 
of £. For each of the models, V(£) — > as £ — > oo. We now consider a photon with 
finite and positive £ which approaches the source from £ = oo, reaches a minimum 
distance from the source at £ m i„ and then returns to £ = oo. In a general spacetime 
the deflection angle is given by 

A = 2^°° ^-n (67) 

4-2.1. AQUAL From the definition of L, as before, we may obtain an expression for 
d^/dip, which turns out to be: 

rf£ _ e 2 *Ly/2(£ - V) 
d^ ~ 2V 
The deflection angle is then given by 

/OO T 
1 ^=^= dr-Tr (69) 

The radius r m i n of closest approach is simply found from the condition that £ = at 
that point, and is found to be r m i n — y/L 2 /(2£) Evaluation of the above integral yields 
A = 0. This is to be expected. The asymptotic AQUAL metric is conformally flat, 
and it is well known that the lensing of photons is insensitive to conformal rescalings 
of the metric. It is this property that leads to substantial conflict with data as the 
deflection of light around galaxies is known to strongly indicate the presence of mass 
discrepancies. The conditions under which a hypothetical purely metric theory of 
gravity underlying MOND would lead to this problem with the asymptotic metric 
have been explored in detail in |36) . 

4-2.2. TeVeS We may proceed as before. The equation for d^/dip again takes the 
form (|68p. However, unlike the AQUAL case, £ is identical to r and the functional 
forms of e 2 * differ. 

The deflection angle is then 

L I r \ 1 



A = 2 — — dr - 7T (70) 




r 2 \ ro J 
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Where we have found r min using the previously described approach, yielding: r min — 
i 

r (2tq£/L 2 ) 2 < 2c -!) . For evaluation of the integral it is convenient to introduce the 

new variable y = L/(roV / 2£)(r/ro) 2c_1 . In terms of this the integral becomes: 

2 /■"-•• 1 J 
A = / — dy — n 

If c < 1/2 then y ma x = 0, y m in = 1, m which case we have that A = 2c7r/ (1 — 2c). 
Therefore the deflection angle is nonzero and constant across different values of r m i n . 
This can be compared to the velocity of circular orbits and precession rate of almost 
circular orbits which also are independent of the distance from the source characteristic 
to the problem. If c = 1/2 then V(^) becomes position independent and so there is 
nothing to prevent the photon falling towards smaller radii than the modified gravity 
regime covers. If 1 > c > 1/2 then via our classical mechanics analogy V{£) no longer 
represents a repulsive force and so there is simply nothing to stop a photon impinging 
from £ = oo from heading to smaller values of £ than are covered by the asymptotic 
metric- i.e. the photon inevitably leaves the modified gravity regime. 

5. Conclusions 

In this article we have considered the properties of timelike and null geodesies for a set 
of metrics. Though these metrics follow from specific models, we anticipate that they 
display features which are not merely confined to those models. It is conceivable that 
any of the metrics may follow from a conceptual framework that does not even involve 
a modified Poisson equation in the transition from General Relativity to modified 
gravity regime. 

By considering a fully relativistic treatment of possible asymptotic geometries for 
a theory of modified Newtonian dynamics, we have found counterparts to the more 
familiar 'flat rotation curve' property of the geometry- namely that the deflection 
angle for null geodesies and the precession rate of near circular orbits may additionally 
become independent of distance from the source. Furthermore, we have in the case 
of TeVeS that the 'scale invariance' inherent to the asymptotic geometry leads to a 
sensitivity to mass that does not exist in General Relativity, specifically that exotic 
behaviour emerges whenever c = VGatMoo > 1, for any value r ^> r , namely that 
the spacetime exchibits a physical singularity as r — > oo for c > 1. We interprete this 
result as saying that physical spacetimes in TeVeS exists only if < \/GNMa n < 1. 

A common property of all the geometries considered here is that they all forbid 
unbound orbits for massive particles (the escape velocity is the speed of light) and all 
of them will lead to a total diminishment of a test photon's frequency as it reaches 
spacelike infinity. 

Additionally, we have seen in the form of the GEA model that a consistent 
relativistic model of MOND exists which agrees with nonrclativistic phenomenology 
but which differs markedly in its behaviour from the other models as one approaches 
a radius r a . For a source of mass M, the ratio of this radius to the radius at which the 
nonrelativistic modified gravity becomes dominant is of order e 3 /VGMa _ j n ^ ne mos ^ 
general terms then it seems difficult then to anticipate the geometry of an underlying 
theory of MOND as r -» oo. 
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